Are Quantum-Cascade Lasers really quantum? 
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The first global simulation of semiconductor-based quantum-cascade lasers is presented; Our fully 
three-dimensional approach allows to study in a purely microscopic way — without resorting to phe- 
nomenological parameters — the current-voltage characteristics of state-of-the-art unipolar nanos- 
tructures. Based on the proposed theoretical scheme, we are able to give a definite answer to the 
long-standing controversial question: is charge transport in quantum- cascade lasers mainly coher- 
ent or incoherent? Our analysis clearly shows that a proper inclusion of carrier-phonon as well as 
carrier-carrier scattering within a semiclassical framework gives excellent agreement with experi- 
mental results. 
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Since the seminal paper of Esaki and Tsu [Q, 
semiconductor-based nanometric heterostructures have 
been the subject of an impressive theoretical and exper- 
imental activity, due to their high potential impact in 
both fundamental and apphed research ||,^ . One of the 
main fields of research focuses on exploiting "band-gap 
engineering" , namely the splitting of the bulk conduction 
band into several subbands, to generate and detect elec- 
tromagnetic radiation in the infrared spectral region, as 
already envisioned by Kazarinov and Suris Q. 

Unipolar coherent-light sources like quantum-cascade 
lasers (QCLs) H , are usually modelled in terms of three- 
levels systems |6|]. Their theoretical description is thus 
often grounded on purely macroscopic models: the only 
relevant physical quantities are the various carrier con- 
centrations within the different energy levels, whose time 
evolution is dictated by proper sets of coupled rate equa- 
tions. This simplified three-level picture does not take 
into account the existence of transverse or "in-plane" de- 
grees of freedom and thus provides no information on the 
microscopic processes governing carrier dynamics within 
the three-dimensional (3D) multi-subband structure. As 
pointed out in Q , such a macroscopic modeling can only 
operate as an a posteriori fitting procedure: The best-fit 
phenomenological parameters obtained through a com- 
parison with experiments do strongly depend on details 
of the 3D hot-carrier distribution and therefore on the 
device operation conditions. 

QCLs are complex devices, whose core is a multi- 
quantum-well (MQW) structure made up of repeated 
stages of active regions sandwiched between electron- 
injecting and collecting regions. When a proper bias 
is applied, an electron cascade along the subsequent 
quantized-level energy staircase takes place. For a de- 
tailed understanding of the basic physical processes in- 
volved, two main issues need to be addressed: (i) the na- 
ture of the hot-carrier relaxation within the device active 
region; (ii) the nature — coherent versus incoherent — 



of the physical mechanisms governing charge transport 
through injector/active-region/collector interfaces. 

Point (i) has been recently addressed in |Q, where the 
macroscopic treatment of hot-carrier relaxation within 
the device active region previously mentioned has been 
compared to a fully kinetic description, based on the fol- 
lowing set of equations : 



^/kiy — [.9kiy — Tkiy/kiyli/i +^[-Pki/,k'i/'/k'i/' — -PkV,kiy/ki/] 
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Here, the first two terms describe — still on a partially 
phenomenological level — injection/loss (i/1) of carriers 
with wavevector k in subband v, while the last ones 
describe intra- as well as inter-subband in- and out- 
scattering processes (ki/ — s- kV) within the device ac- 
tive region only. The Boltzmann-like structure of (0) 
allows for a partially stochastic solution; to this aim, the 
commonly used Monte Carlo (MC) method |^ has been 
adopted. The simulated experiments in |Q show that, 
as far as the active region is concerned, the quantum- 
cascade process is mainly determined by carrier-optical 
phonon emission; Moreover, the values of effective inter- 
level rates strongly depend on the shape of the in-plane 
charge distributions as well as on Pauli-blocking effects. 

However, the microscopic analysis in [0, being lim- 
ited to the device active region only, does not al- 
low to answer point (ii): which are the basic phys- 
ical mechanisms governing charge transport through 
injcctor/active-region/coUcctor interfaces? This issue 
is intimately related to the long-standing controversial 
question ||^: is charge transport in quantum- cascade 
lasers mainly coherent or incoherent? 

To provide a definite answer to this fundamental ques- 
tion, in this Letter we present the first global simu- 
lation — injector plus active region plus collector — of 
semiconductor-based QCL structures. To this end, the 
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partially phenomenologic model in (|]) has to be replaced 
by a fully microscopic description of the whole MQW 
core structure. This is not a trivial task. It requires: (i) 
a proper description of all the 3D electron states within 
the active as well as injector/collector regions; (ii) a fully 
microscopic simulation scheme to "close the circuit" . 

As a starting point, we evaluate — within a self- 
consistent Schrodinger-Poisson approach in an envelope- 
function framework — the set of 3D single-particle 
electron states {k:^} corresponding to a single QCL stage, 
e.g., active region plus collector/injector in the presence 
of the appHed bias (see Fig. |^). Given such carrier states, 
we consider the ideal MQW structure obtained as infinite 
repetition of this QCL periodicity region (see Fig. |^). 
Within such extended scheme, the time evolution of the 
carrier distribution function /ka is governed by the fol- 
lowing Boltzmann-like equation: 



dt 



(2) 



Here, a = (A, f) denotes the generic electron state in our 
MQW structure, i.e., the ly-th state of the A-th stage, 
while s labels the different scattering mechanisms con- 
sidered, e.g., carrier-phonon, carrier-carrier, etc. 

To "close the circuit" , we impose periodic boundary 
conditions limiting the inter-stage (A' 7^ A) scattering to 
just nearest-neighbor coupling (A' = A ± 1). In view of 
the translational symmetry, we are allowed to simulate 
carrier transport over the central — i.e., A = — stage 
only. Every time a carrier in state v undergoes an inter- 
stage scattering process (i.e., 0, — > ±1, i/', the electron 
is properly reinjected into the central region (0, v 0, v') 
and the corresponding electron charge ±e will contribute 
to the current through the device |pl]| . The proposed 
global-simulation scheme allows to study in a purely mi- 
croscopic way — without resorting to phenomenological 
i/1 parameters — the current-voltage characteristics of 
statc-of-thc-art unipolar nanostructures. 

At this point a few comments are in order. Contrary to 
Eq. (|l|), the Boltzmann-like equation in (^) — describing 
an infinite system — does not contain explicit i/1 terms, 
which are automatically accounted for within the pro- 
posed re-injection scheme. Moreover, the electric field 
does not enter our transport equation via the usual drift 
term 10 1; the corresponding bias is directly taken into 
account in evaluating the single-particle states {k;^} [ pT| . 

It is well-known that the application of the present 
semiclassical transport theory to MQW structures at 
high fields becomes questionable. Indeed, according to 
(^, in the absence of scattering processes (Pka,k'Q' — 0) 
we have no current. This paradoxical result is due to 
the neglect of phase coherence between our localized 
states. A proper treatment of coherent transport — 
including resonant-tunneling effects — would require a 
density-matrix formalism 02]. However, in the presence 



of a strong scattering dynamics current is mainly deter- 
mined by incoherent hopping processes between localized 
states; moreover, such phase-breaking processes lead to a 
very fast decay of off-diagonal density-matrix elements. 
In this regime coherent-transport effects play a minor 
role and the semiclassical limit is recovered. It is then 
clear that the nature — coherent versus incoherent — of 
charge transport in QCLs is primarily dictated by the 
role played by incoherent scattering processes. 

In order to discriminate between the two different 
regimes, we have applied the above global-simulation 
scheme to state-of-the-art QCL structures. In particular, 
as prototypical device, we have considered the diagonal- 
configuration QCL in schematically depicted in 
Fig. |l|, in which the proposed simulation scheme is also 
sketched. Here, the energy levels and probability densi- 
ties of various electron states within the simulated stage 
(A = 0) are also plotted: They correspond to the device 
active region (full line, — 1,2,3 according to the stan- 
dard notation) as well as to the collector region (dashed 
line, = A,B,C,D,E). The nearest stages A = ±1 are 
partially shown for clarity as well. 

In order to properly model phase-breaking hopping 
processes, in addition to carrier-optical phonon scatter- 
ing, all various intra- as well as intersubband carrier- 
carrier interaction mechanisms have been considered. To 
this end, the well-established MC method for the simula- 
tion of intercarrier scattering in quantum- well systems |^ 
has been extended to our MQW biased structure. Within 
such approach — based on the Born approximation — the 
effect of nonequilibrium multisubband screening of the 
bare two-body Coulomb potential is taken into account 
within the standard random-phase approximation [Q. 
Other scattering mechanisms, not included in the simu- 
lation, are expected to play a very minor role j^: The 
interaction with acoustic phonons — which in the ab- 
sence of carrier-carrier scattering would be essential to 
get ergodicity — does not affect energy relaxation signifi- 
cantly due to its quasi-elastic nature and small coupling 
constant; Carrier- impurity scattering is also expected to 
have very little impact on vertical-transport phenomena. 

As a first step, we have investigated the relative 
weigth of the carrier-carrier and carrier-phonon compet- 
ing energy-relaxation channels. The time evolution of 
the carrier population in the various subbands as well 
as of the total current density in the presence of carrier- 
phonon scattering are depicted in Fig. 2. Parts (a) and 
(b) refer, respectively, to simulated experiments with- 
out and with two-body carrier-carrier scattering. In our 
"charge conservation" scheme, we start the simulation 
assuming the total number of carriers to be equally dis- 
tributed among the different subbands; then the electron 
distribution functions evolve according to Eq. (^ and 
a steady-state condition is eventually reached, leading 
to the desired 3 — > 2 population inversion. As we can 
see, the inclusion of intercarrier scattering has signifi- 
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cant effects; It strongly increases inter-subband carrier 
redistribution, thus reducing the electron accumulation 
in the lowest energy level A and optimizing the coupling 
between active region and injector/collector (the popula- 
tions of subbands 3 and B get equal). This effect comes 
out to be crucial in determining the electron flux through 
the MQW structure. Indeed, Fig. 2(c) shows the simu- 
lated current density across the core region of our QCL 
device, obtained with (full line) and without (dashed line) 
the inclusion of carrier-carrier scattering. These simu- 
lations have been performed at the threshold operating 
parameters estimated in |Q . The current density in the 
presence of both electron-phonon and electron-electron 
scattering mechanisms is about 4 kA/cm^: This value is 
in good agreement [p^ with the experimental results in 

The results plotted in fig. g clearly demonstrate that 
the electron-phonon interaction alone is not able to effi- 
ciently couple the injector subbands to the active region 
ones: While carrier-phonon relaxation well describes the 
electronic quantum cascade within the bare active re- 
gion [Q, carrier-carrier scattering plays an essential role 
in determining charge transport through the full core 
region. This can be ascribed to two typical features 
of carrier-carrier interaction — compared to the case of 
carrier-phonon — : (i) this is a long-range two-body in- 
teraction mechanism, which also couples non-overlapping 
single-particle states [see Fig. 1]; (ii) the correspond- 
ing scattering process at relatively low carrier density is 
quasielastic, thus coupling nearly resonant energy levels, 
like states 3 and B. 

This result is extremely important: it tells us that 
the presence of a coherent-tunneling dynamics — often 
invoked to account for nearly resonant inter-level pro- 
cesses between active region and injector/collector — is 
not necessary: The charge transport in QCL structures 
can be explained quantitatively in terms of incoherent 
intercarrier-assisted tunneling processes, and thus within 
a purely semiclassical transport theory. This allows us to 
answer the long-standing controversial question on the 
nature — coherent versus incoherent — of charge trans- 
port in QCLs previously mentioned: for the typical struc- 
tures considered, energy-relaxation and dephasing pro- 
cesses are so strong to destroy any phase-coherence ef- 
fect; as a result, the usual semiclassical description is in 
excellent agreement with experimental results. 

Since in the proposed MC simulation scheme, the cur- 
rent density across the whole structure is an output 
of the simulation, the current-versus-voltage behavior 
of our semiconductor quantum device can be obtained 
within our fully microscopic description, without resort- 
ing to any phenomenological external parameter. Figure 
3 presents the simulated voltage-current characteristics 
for the MQW structure of Fig. As we can see, the cur- 
rent density increases with the applied field, due to the 
increased coupling between subband 3 and various lower- 



energy injector states, in which carriers tend to relax and 
accumulate. The agreement between our theoretical re- 
sults and the measured dependence shown in Ref. |l^] 
is extremely good. Discrepancies between simulated and 
measured values are well within the measurement uncer- 
tainties and the complexity of the real device structure. 

In summary, we have presented the first global simu- 
lation of semiconductor-based QCLs. The proposed the- 
oretical scheme allows for a fully microscopic evaluation 
of the current-voltage characteristics of prototypical QCL 
structures. The excellent agreement between our semi- 
classical simulations and recent experimental results con- 
firms the incoherent nature of charge transport in such 
unipolar quantum devices. 

This work has been supported in part by the European 
Commission through the Brite Euram UNISEL project 
(Contract No. BE97-4072). 
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FIC. 1. Schematic representation of the conduc- 

tion-band profile along the growth direction for the diago- 
nal-configuration QCL structure of Ref. [^. The MQW is 
biased by an electric field of 48 kV/cm. The levels = 1, 2, 3 
and 1/ = A,B,C,D,E in the active and collector regions (full 
and dashed lines, respectively) of the simulated stage (A — 0) 
are also plotted together with the corresponding probability 
densities. The replica of level 3 in the following stage A = -1-1 
is shown for clarity. 

FIG. 2. Time evolution of simulated carrier densities in 
the various subbands of the MQW structure of Fig. ^ (full 
lines: f — A,B,C,D,E; dotted lines: v — 1,2,3), without 
(a) and with (b) intercarrier scattering, (c): current density 
across the whole structure in presence (full line) and absence 
(dotted line) of carrier-carrier interaction. An electric field of 
48 kV/cm has been applied. 

FIG. 3. Applied-field vs current-density characteristic at 
77K of the MQW structure of Fig. 0. The dashed hue is a 
guide to the eye. 
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R. C. lotti and F. Rossi, Figure 1 




R. C. lotti and F. Rossi, Figure 2(a) 




R. C. lotti and F. Rossi, Figure 2(b) 




R. C. lotti and F. Rossi, Figure 2(c) 
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R. C. lotti and F. Rossi, Figure 3 



